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Abstract 

Differences among calculations of the atmospheric neutrino beam 
can be traced in large part to differences in the representation of pion 
production by protons interacting with nuclei in the atmosphere. In 
this paper we review the existing data with the goal of determining 
the regions of phase space in which new measurements could help to 
improve the input to the calculations. 

1 Introduction 

The observed up-down asymmetry of the flux of muon-hke events induced 
by interactions of atmospheric neutrinos has been interpreted as evidence for 
neutrino oscillations The case for neutrino oscillations is especially strong 
for the multi-GeV events at Super-Kamiokande because the energies are high 
enough so that the charged leptons follow the neutrino direction closely and 
expose the oscillation effects as a function of the neutrino pathlength. On 
the other hand, comparisons to calculations of the atmospheric neutrino flux 
over as broad an energy range as possible are essential for making a detailed 
interpretation of the data and generally for giving confidence in the basic 
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result. In addition, calculations play a crucial role in comparison of measure- 
ments made at different locations and with different techniques. In particular, 
the Soudan experiment 0|, and 1MB previously, detect neutrino interac- 
tions in a range similar to Super-Kamiokande (and Kamiokande [H), but the 
geomagnetic environments are quite different, which changes the neutrino 
angular distribution expected in the absence of oscillations. Measurements 
of neutrino- induced upward muons, as also at MACRO 0, are also sensitive 
to oscillations, but in a different energy regime (see Fig. |l]). 
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Figure 1: Distributions of neutrino energies that give rise to four classes of 
events. Sub-GeV and multi-GeV refer to the two classes of contained events 
at Superkamiokande. The response for contained events is averaged over all 
angles, while for stopping and throughgoing muons it is only given for vertical 
upward going neutrino induced muons. 

There are several independent calculations of the flux of atmospheric 
neutrinos [0-0]. In all cases but one |^, the calculations start with the 
primary cosmic-ray spectrum as measured (protons and alpha-particles give 
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the dominant contributions) and simulate the resulting cosmic-ray cascade 



in the atmosphere, calculating the fluxes of v^. 



and /i . The 



calculated muon fluxes can be compared with measurements of muons at 
high altitude |[T^-|T6| as a check of the calculation. Perkins ||^ starts from 
the measured muon flux and uses the kinematic relationship between muons 
and neutrinos to derive the neutrino flux. A limitation of the muon fluxes is 
that the individual measurements generally have low statistics and depend 
more sensitively than the neutrinos on the details of geomagnetic and three- 
dimensional effects. 

The most recent calculations 



10,11 



constitute a very signiflcant tech- 
nical advance in that they are three-dimensional. In previous calculations 
secondary particles had been assumed to follow the direction of the primary 
cosmic-ray that generated them. A major conclusion of Ref. |0 is that for 
practical purposes the one- dimensional calculations give adequate results. 
This is because Fermi-momentum of the target nucleons of the neutrino in- 
teractions, compounded by limited angular resolution, smear out angular 
features visible in fluxes of ~GeV neutrinos. A corollary of this conclusion 
is that simpler one-dimensional calculations can be used for comparison of 
calculations in order to trace sources of differences among the calculated 
neutrino fluxes. 



Some time ago the existing one-dimensional calculations were compared JT? 
The dominant sources of difference were found to be the assumed primary 
spectrum and the representation of pion production in collisions of protons 
with light nuclei. Although differences in input were at the level of 20-30% 



- or larger in some cases - the neutrino fluxes of the two calculations |18|,[19 



that have been used extensively for input to the detector simulations differ 
by much less. This is a consequence of compensating uncertainties together 
with the fact that both calculations were constrained to flt the same (ground- 
level) muon data. Recently, several new measurements of the primary spec- 
trum |E0|-E2l conflrm the lower |E3| of two older measurements P3|,|2^, leav- 



ing differences in representation of pion production as the major source of 
uncertainty for GeV neutrino fluxes. 

A striking example is the difference between the calculations of Battistoni 
]TU[ and the Bartol fluxes H. This comparison, which was made P3 



et al. 



T0[| and the same assumed primary 

for 



using the one- dimensional version of Ref 

spectrum as Ref. 0, gives neutrino fluxes 20-30% lower than Ref. 

< 10 GeV. A detailed comparison with the new calculation of Ref. |]I0| is 
currently in progress |26]. 
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2 Inclusive pion production around 20 GeV 



To see what range of interaction energies is most important for the vari- 
ous classes of events shown in Fig. |l|, one needs to look at the distribu- 
tions of primary energies that produce the events. Roughly speaking, for the 
steep cosmic ray spectrum the most relevant primary energies are an order 
of magnitude higher than the neutrino energy of interest. For the upward 
neutrino-induced muons, for example, the important range of interaction en- 
ergies extends up to several TeV. In this energy range, the uncertainties in 
the primary spectrum are still relatively large, and there are also significant 
uncertainties from the amount and momentum distribution of kaon produc- 
tion 1 28], which will be addressed in [2^. Here we concentrate on the lower 



energy events (e.g. the sub-GeV events at Super-Kamiokande) for which 
most of the contribution comes from interactions of primary cosmic-ray nu- 
cleons with energies between 10 and 100 GeV. A similar range of energies 
is responsible for the Soudan events and for contained events in 1MB and 
Kamiokande. 

We show the distribution of primary energies that gives rise to the sub- 
GeV events at Super-Kamiokande in Fig. |^ [^. We show separately the 
response for downward (C) and upward (B) events, as well as the response 
that would apply if there were no geomagnetic field (A). It is interesting to 
note that the geomagnetic cutoff leads to an observable site-dependence of the 
ratio of downward to upward moving events. In the absence of oscillations, 
due to the higher local geomagnetic cutoff at Kamioka in contrast to the very 
low vertical cutoff in the northern U.S., the ratio of downward to upward 
going leptons is respectively less and greater than one. 

One of the most extensive data sets in the energy range responsible for 



sub-GeV events is the work of Eichten et al. p9[. Figure 3 shows the data 
from this experiment for production of vr"*" from interactions of 24 GeV/c 
protons on beryllium. The solid lines are fits of the form 

'^^ dpd^pT ^ "^^^ ^ {-Kp) ^t} , (1) 
where the transverse mass is defined as 



mT = \/pl + ml, (2) 

and a{p), b{p) are free parameters at each value of pion momentum p. This 
form gives good fits for all momenta except the lowest two values. The dashed 
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Figure 2: Response for sub-GeV muon-like events in Super-K to the energy 
of the primary cosmic ray nucleons. A: no geomagnetic cutoff; B: events from 
lower hemisphere (upward going leptons); C: events from upper hemisphere 
(downward going leptons). Each pair of curves shows the range of the signal 
between minimum (solid) and maximum (dotted) of solar activity. 

lines show fits in which the pion mass entering Eq. (Q) is also treated as a free 
parameter. For negative pions, the form ([l|) gives good fits for all measured 
momenta. 

What is important for the overall normalization of the atmospheric neu- 
trino flux is the total yield of pions integrated over all phase space. This in- 
volves both interpolation and extrapolation. Figure ^ illustrates the fraction 
of the pion production cross section that requires interpolation only for the 
data of Ref . |^ . The curve represents the fraction of secondary pions falling 
into the acceptance range of the experiment as function of a;iab = p/Pheam, 
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Figure 3: Data on inclusive vr"*" production in proton-beryllium collisions at 
24 GeV The solid lines correspond to fits using the pion mass in and 
the dashed curves show the result if the mass is treated as free parameter. 



using the predictions of the TARGET Monte Carlo [pO|] . 

Thus, for example, for 7r+ at p = 6 GeV (xiab ~ 0.25) approximately 75% 
of the integral over transverse momentum does not require extrapolation 
into unmeasured regions of transverse momentum. Most important is the 
fact that in this experiment there are no data at all for p < 4 GeV. In this 
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Figure 4: Fraction of pion production cross section covered by the data of 
Eichten et al. (p-Be at 24 GeV Allaby et al. (p-Be at 19.2 GeV 

and Abbott et al. (p-Be at 14.6 GeV ||33|). For Eichten et al. the upper 
and lower curves show the coverage for vr^ and vr" respectively. The other 
experiments have acceptances which are essentially identical for negatively 
and positively charged pions. 



region, therefore, the representation of pion production depends on how well 
the model extrapolates in both longitudinal and transverse momentum. 

Other experiments |pTl-33| cover various regions of phase space at various 
nearby beam energies (see Fig. |). The Lundy et al. data set |^ (p-Be 
at 12.5 GeV) is not shown since it corresponds to pion production spectra 
which are not consistent with the other measurements. 

Only the data of Abbott et al. |^3| reach into the xiab region of about 0.1 
to 0.2, which is an important contributor to the atmospheric lepton fluxes. 
Since these measurements are published as function of particle rapidity y = 
log[(£' + pl)/{E — Pl)] and transverse mass m^, the conversion of the data 
into dN/dx\ah requires a two-dimensional inter- and extrapolation resulting 
in substantial uncertainties. 
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In Fig. 1^ we show the inclusive pion spectra obtained by integrating ana- 
lytic parametrizations fitted to the data of Eichten et al. |29|, Allaby et al. 
and Abbott et al. I^B] . The sensitivity of the results to different parametriza- 
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Figure 5: Inclusive pion production spectra obtained from Eichten et al. |2 
(solid circles), Allaby et al. |^ (triangles), and Abbott et al. (diamonds) 



tions used for extrapolation is indicated by showing two points at the same 
X- value for a given data set. In addition to the typical systematic uncertain- 
ties of these experiments of about 15%, different methods of extrapolation 
result in pion spectra which are different by up to 25%. 

Finally, there are many proton-beryllium and proton-carbon particle pro- 
duction measurements published which are restricted to one or two angular 



bins of the secondaries (for example. Baker et al. ||3^, 10, 20, 30 GeV; Dekkers 
et al. 19, 23 GeV). Due to the lack of theoretical understanding of soft 
hadronic particle production these data sets cannot be used for inter- and 
extrapolation to obtain secondary pion spectra. 



3 Pion production other energies 



There are several measurements at lower energy, which cover some relevant 
regions of phase space [p^-|47[|. In most of the cases only vr" distributions 
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have been measured and some of the papers report on mean multiphcities 
or multiphcity distributions only. Also, many data sets refer to thick tar- 
gets and are not directly usable for cascade calculations. In the following 
we discuss only some of the published data sets where final state particles 
have been identified and a large fraction of the phase space was covered by 
the experiment. In particular, we do not consider the numerous emulsion 
measurements. 

At low energy, the data of |]39|-p| are interesting because of the full 
coverage of the phase space. For example, in Fig. ^ we show vr" produc- 
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Figure 6: Inclusive production spectrum of negative pions in proton-carbon 
collisions [B^,^ . 



tion spectra for proton-carbon collisions at 4.2 and 10 GeV. The mean vr 
multiphcities are respectively 0.33 ± 0.02 and 1.00 ± 0.03 In addi- 

tion, in Refs. PD| , |^ data on a-carbon and carbon-carbon collisions with 4.2 
GeV/nucleon are given. 

The most relevant data sets at higher energy known to the authors are in 



Refs. 



Ref. ||4^ gives the ratio of particle production yields in proton- 
beryllium and proton-aluminum collisions at 67 GeV. Such studies are im- 
portant for the extrapolation of p-Be data to p-air collisions. Whereas the 



-'^In Ref. [Q the 4.2 GeV data are also given as double-differential distributions 
Ed^a/dp^ in 10 degree bins of the azimuthal angle. 
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are 



measurements published in [JS| refer to thin targets, the data of Ref. |^ 
for thick targets and both thick and thin targets are considered in Ref. |49 



The measurements of p-Be colhsions at 450 GeV of the NA56/SPY collabo- 
ration cover the important transverse momentum range of — 600 MeV, 
however cross sections obtained by extrapolation to a thin Be target are given 
only for the very forward direction. 



4 Conclusion 

After the new measurements of the primary cosmic ray spectrum, the un- 
certainties in the calculated sub-GeV neutrino fluxes are dominated by the 
limited information and understanding of hadronic interactions in the energy 
range from about 5 to 100 GeV. 

Most pressing is the need for a single experiment that covers pion produc- 
tion at several beam energies in the peak region of the sub-GeV events in Fig. 
2, namely, around 20 GeV/c. Targets should be as close as possible to the 
constituents of the atmosphere, especially nitrogen. If this is not possible, a 
series of targets with mass number spanning A=14 should be done. 

At low energy, several data sets cover vr" production. Data on 7r+ pro- 
duction is very sparse and many experiments refer to thick targets or suffer 
from low statistics. Thus, a consistent set of new measurements extending 
to lower energy would also be important, especially for interpreting measure- 
ments made at low geomagnetic cutoff. 

Use of a beam of helium nuclei would also be of interest. In the en- 
ergy range from 10 to 100 GeV/nucleon, approximately 80% of the primary 
cosmic-ray nucleons are free protons and 15% are bound in alpha particles. 
Only 5% or less are in heavier nuclei. 
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